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LS 5039 is a TeV gamma-ray binary with extended radio emission. It consists of a compact ob- 
ject in the mildly eccentric (e=0.35), 3.9-day orbit around a massive O star. The nature of the 
compact object is not yet established. In this paper, assuming that the compact object is a black 
hole, we study the accretion of 0-star wind by the black hole, by performing three-dimensional 
Smoothed Particle Hydrodynamics (SPH) simulations. In order to roughly emulate the effect of 
the stellar radiation effectively canceling the stellar gravity, we assume that the O star's gravity 
does not exert on the wind. The wind particles are ejected with half the observed terminal veloc- 
ity in a narrow range of azimuthal and vertical angles toward the black hole, in order to emulate 
the wind significantly slower than the terminal speed, and optimize the resolution and computa- 
tional efficiency of simulations. We find that the mass-accretion rate closely follows the classical 
Bondi-Hoyle-Littleton accretion rate, which is of the order of lO^^gs"^ around periastron. The 
accretion rate at this level would provide jets enough power to produce the gamma-rays detected 
by HESS. Since the accretion peak occurs near the periastron passage, we need a strong gamma- 
ray absorption around periastron in order for the microquasar scenario to be consistent with the 
observed orbital modulation of the TeV gamma-ray flux. 
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1. Introduction 



The high mass X-ray binary LS 5039 consists of an 06.5V star and a compact object in the 
^ 3.9-day orbit with a moderate eccentricity (e ^ 0.35). It is one of the three binaries from which 
persistent TeV y-rays have been detected [jl|]. The VHE y-ray luminosity from LS 5039 modulates 
between 4 — 10 x 10^^ ergs~^ The maximum flux takes place around phase 0.8-0.9 ^ ||, 0], 
where phase corresponds to the periastron passage. The system shows a persistent radio emission, 
the morphology of which is consistent with mildly relativistic jets [14, 15]. The nature of the 
compact object has not been established. proposed that it is a black hole, which would place 
LS 5039 in the microquasar class, whereas some argue that the system is a colliding-wind binary 



with a non-accreting pulsar (e.g., Qq, 11, 17]). To test the microquasar scenario for LS 5039, it is 



important to study whether the black hole can get high enough mass-accreion rate to power jets 
to produce the y-rays detected by HESS. In this paper, we briefly report on the result from three 
dimensional, numerical simulations of wind accretion by a black hole in LS 5039. 



2. Numerical model 



Simulations are performed with a 3-D Smoothed Particle Hydrodynamics (SPH) code. The 
code is based on a version originally developed by Benz ^ Q], then by Bate and his collaborators 

and recently by [ |T3[ ] and [|T^. It uses the variable smoothing length and the SPH equations 
with the standard cubic-spline kernel are integrated with individual time step for each particle. The 
artificial viscosity is assumed to have the standard values of aspH = 1 and jSsph = 2. In our code, 
the O-star wind is modeled by an ensemble of isothermal gas particles of negligible mass, while 
the compact object, which we assume to be a black hole, by a sink particle with the corresponding 
mass and a radius r^cc taken to be much smaller than the classical Bondi-Hoyle-Littleton (BHL) 
accretion radius, tbrl. at periastron. If gas particles fall within this radius, they are removed 
from the simulation. On the other hand, we model the O star as a particle with the corresponding 
mass, whose gravity exerts only on the binary companion, not on the wind particles. We take this 
assumption to roughly emulate the effect of the stellar radiation effectively canceling the stellar 
gravity. Hence, in our simulation, the gravity by the compact object is the only external force 
exerting on the wind particles. The wind particles have the initial velocity of half the observed 
terminal velocity, given that the accretion onto the black hole occurs in the region where the wind 
speed is still significantly slower than the terminal velocity. We set the binary orbit in the x-y plane 
and the major axis of the orbit along the x-axis, where the apastron is in the +x-direction. The 
stellar, wind, and orbital parameters adopted for LS 5039 simulations are summarized in Table [T[ 

To optimize the resolution and computational efficiency of our simulations, the wind particles 
are ejected only in a narrow range of azimuthal and vertical angles toward the black hole ^ We 
have confirmed that this method provides a quite similar flow structure to that of a corresponding 

^In the poster presented at the Workshop, we had broken the stellar wind and accretion portions of the computation 
into two separate, but linked parts, in which the first focused on the stellar wind expansion and the rate of mass that 
enters a spherical region around the balck hole, and the second focused the mass accretion onto the black hole, using the 
information of the first part of the simulation. However, after the Workshop, it turned out that the second part had a bug 
that seriously affected the final result. Therefore, we have revised our code and this time rerun simulations of the whole 
system, using a much smaller accretion radius. The result reported in this paper is based on these revised simulations. 
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Table 1: Model parameters for LS 5039 



Primary 



Secondary 



Spectral type 
Mass (M0) 
Radius 

Effective temperature Tgff (K) 
Initial wind velocity yw,ini (kms~^) 
Mass loss rate (M0yr"^) 
Orbital period Porb (d) 
Orbital eccentricity e 
Semi-major axis a (cm) 



06.5V Black Hole 

22.9^ 3.T 
9.3/^0^= 0.31a) 2.5x 
39,000^ 
1,200 
5 X 10-^ ^ 

3.9060^ 
0.35^ 
2.17 X 10^2 



The upper limit of the low X-ray state and the lower limit of the high state [0]. 



spherically symmetric wind. Moreover, the former has a much higher spatial resolution than the 
latter, as expected. We can see this in Fig. |l| comparing the density distribution in the orbital plane 
at periastron from such a simulation with that from a simulation of an axisymmetric wind ejected 
within a constant opening angle of 20°. Both simulations use similar number of wind particles 
A^SPH- In the former simulation, the accretion radius of the black hole r^^c = 2.5 x 10~^, while 
^acc = 10~^ for the latter simulation because of the lower spatial resolution. 



D 




D 





Figure 1: Wind density in the orbital plane at periastron. The left panel shows the logarithmic density 
from a simulation where wind particles are ejected toward the black hole (a small bright dot on the left of 
the O star) in the restricted range of 30° in the azimuthal direction and 20° in the polar direction, whereas 
the right panel shows that from a simulation of an axisymmetric wind with the same opening angle of 
20°. Annotated in each panel are the phase, the time from the periastron passage, and the number of wind 
particles, race = 2.5 x 10~^ in the left panel, while race = 10~^ in the right panel because of the lower spatial 
resolution. 



3. Results 

In Figs. 1^ and ||, we present the density distribution of the wind in the binary orbital plane at 
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several different phases (0 = corresponds to periastron). Figure ^ covers the whole simulation 
volume (within 1 .5a from the center of mass of the binary), while Fig. || shows the detailed structure 
of the wind around the black hole. In both figures, the color-scale plots show the logarithmic density 
of the wind in the orbital plane, while arrows superposed on the density plot in Fig. || denote the 
velocity vectors of the flow. In Fig. ^ the dashed circle denotes the classical BHL accretion radius, 
^BHL, for a supersonic flow, which is given by 



^BHL 



2GMx 



(3.1) 



•^rel 



where Mx is the mass of the black hole and Vrei is the relative velocity between the wind and the 

-1 



black hole (e.g., [Q]). Here we have calculated tbhl by taking Vw,ini = l,200kms as the wind 
velocity. The line fragment on the dashed circle in Fig. || shows the direction of the O star. 
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Figure 2: Snapshots at six different phases. The color-scale plot shows the logarithmic density in the orbital 
plane in cgs units. The large filled circle near the center is the O star, while a tiny bright dot (with the radius 
of 2.5 X I0~^a) at the tip of bow shocks is the black hole. 

Figure |2| shows that strong bow shocks form behind the black hole and that its shape is distorted 
around periastron by the rapid orbital motion of the black hole. Note that, despite the spiral shape 
of the wind region, which is a direct consequence of our wind ejection method, the flow is radial 
except in a region near bow shocks. 

From Fig. |3[ we can see that the accretion is basically of classical BHL type. The wind passing 
through a region within ^ tbrl from the black hole is strongly focussed and ultimately accreted. 
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Figure 3: Detailed flow structure around the black hole at the same phases as in Fig. 2. Color- scale plot 
shows the logarithmic density in the orbital plane, while arrows denote the velocity vectors of the flow. 
The dashed circle around the black hole denotes the Bondi-Hoyle-Littleton accretion radius, tbhl, for a 
l,200kms~^ wind. The line fragment on the dashed circle shows the direction of the O star. 



while the flow is almost parallel in an outer region. In order to see this more clearly, we compare 
in Fig. the simulated mass-accretion rate, M, with the BHL accretion rate, Mbhl, given by 

Mb„l = ^, (3^2) 

where d is the binary separation (e.g., [Q]). Here we have calculated Vrei, assuming the constant 
wind velocity of Vw,ini- In this figure, we have averaged the rapidly fluctuating accretion rate over 
the phase interval of 0.01 (histogram in blue). The red line denotes the BHL accretion rate, Mbhl- 
It is immediately observed that the simulated accretion rate closely follows the BHL rate. Thus, the 
classical BHL approximation also basically hold in 3-D simulations as it does in 2-D counterparts. 
A closer look at the figure, however, tells us that the simulated accretion rate is slightly lower 
than the BHL rate. This is due to the effect of wind acceleration by the gravity of black hole; 
the assumption of constant wind velocity results in an slightly underestimated relative velocity, 
which in turn leads to a slightly overestimated BHL accretion rate shown in Fig. ^. The simulated 
accretion rate modulates by a factor of ^ 5. It varies from ^ lO^^gs"^ around periastron to ^ 
2 X lO^^gs"^ around apastron. 
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Figure 4: Mass-accretion rate onto the black hole, which is averaged over the phase interval of 0.01. For 
comparison, the BHL accretion rate for a wind of 1, 200 km s~^ is shown by the red line. 



It is important to note that the mass-accretion rate of ^ lO^^gs"^ would provide jets enough 
power to produce the y-rays detected by HESS. With M ^ lO^^gs"^, the accretion power is ~ 
10^^ erg s~^ Assuming that 10% of this power goes to jets and ^ 1 /3 of the jet power is converted 
to relativistic particles, then we expect a y-ray luminosity of roughly ^ 3 x lO^^ergs"^ which 
would explain the observed y-ray luminosity of 2.7 x lO^^ergs"^ above 100 Me V [p^. 

Now, with the accretion rate as high as ^ lO^^gs"^, why do we see no thermal feature at 
X-rays (e.g., [^)? One possibility is that the accretion flow in LS 5039 is radiatively inefficient. 
For a black hole of 3JMq, the Eddington luminosity is ~ 4.6 x lO^^ergs"^ With the energy 
conversion rate ^ 0.1, the critical mass accretion rate, Merit is ^ 5 x 10^^ gs~^ Then, the accretion 
rate of ^ lO^^gs"^ corresponds to ^ 2 x 10~^Mcrit. This rate is well below the maximum accretion 
rate, 0.05 — O.lMcrit, below which the advection-dominated accretion solutions exist [|T2[]. Thus, 
our model is consistent with ADAF solutions, which could explain the paucity of thermal X-ray 
emission in LS 5039. 



4. Conclusions 

We have carried out 3-D SPH simulations of the wind accretion in LS 5039, assuming that 
the compact object is a 3JMq black hole. Wind particles were ejected at half the observed ter- 
minal velocity to emulate the stellar wind significantly lower than the terminal velocity. In order 
to optimize the resolution and computational efficiency of simulations, we have simulated a wind 
restricted to a narrow range of angles toward the black hole. We have found that the simulated ac- 
cretion rate closely follows, but is slightly lower than, the classical Bondi-Hoyle-Littleton accretion 
rate. The maximum accretion rate of ^ lO^^gs"^ for the stellar mass-loss rate of 5 x lO~^M0yr~^ 
takes place at a phase slightly after the periastron passage. This rate would provide jets in LS 5039 
enough power to produce the y-rays detected by EGRET and HESS. Since the accretion peak occurs 
near periastron, we need a strong y-ray absorption around periastron in order for the microquasar 
scenario to be consistent with the observed peak phase of 0.8-0.9 of the TeV y-ray flux. 
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